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SUMMARY 

The deformation at the core-mantle boundary produced by the 2004 Sumatra earthquake is 
investigated by means of a semi-analytic theoretical model of global coseismic and postseis- 
mic deformation, predicting a millimetric coseismic perturbation over a large portion of the 
core-mantle boundary. Spectral features of such deformations are analysed and discussed. The 
time-dependent postseismic evolution of the elliptical part of the gravity field ( J 2 ) is also com- 
puted for different asthenosphere viscosity models. Our results show that, for asthenospheric 
viscosities smaller than 10 18 Pa s, the postseismic J 2 variation in the next years is expected to 
leave a detectable signal in geodetic observations. 

Key words: core-mantle boundary - deformation - rheology - seismic modelling - spectral 
analysis. 



1 INTRODUCTION 

The devastating megathrust earthquake occurred on December 26th, 2004 off the west coast of 
northern Sumatra was the second-greatest event ever registered, according to current estimates 
which give a moment magnitude M w = 9.3. Such an exceptional event produced measurable ef- 
fects on many geophysical observables. It has been shown (|Park et al. 2005|) that the 2004 Sumatra 
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earthquake was able to excite Earth's free oscillations of exceptionally large amplitude. These 
oscillations, characterized by periods T > 1000 s, remained observable for weeks in broad- 
band seismic data providing information on the size and duration of the event itself. At the 
same time, this event is expected to have produced a jump in the rotational poles secular mo- 
tion (|Chao & Gross 20051 IGross & Chao 20061) . yet current geodetic measurements were not able 
to detect it, probably because it is shadowed by other effects of atmospheric and oceanic origin. 
Associated static surface deformation field data shows (Banerjee et al. 2005 ; iBoschi et al. 2006|) 
that the Sumatra earthquake produced static offsets of the order of 1 mm recorded by continuous 
GPS stations located up to 5000 km away from the epicentre. On the basis of the aforementioned 
points, it is reasonable to expect that the static deformation associated with the event could have 
affected also a large part of the Earth's interior. 

The main aim of this work is to contribute to the characterization of the global effects of Suma- 
tra event on core-mantle boundary (CMB) topography and on the J 2 gravitational field coefficient. 
Using a semi-analytical deformation model, we estimated the amplitude and shape of the CMB 
topography changes as the result of this event both in the purely elastic (t = 0) and fluid (t — > oo) 
limits. In addition, we have performed a spherical harmonic decomposition of the coseismic CMB 
deformation field, investigating its symmetric properties. This analysis has been motivated by the 
possible connection between CMB deformations with axial and equatorial symmetry and core flow 



perturbations (Dumberry & Bloxham 2004). As a result, we have found that the axial and equa- 
torial symmetric component of the CMB deformation has an amplitude of the order of a fraction 
of millimeter. It is worth noting that, although the effect of Sumatra earthquake on the CMB may 
appear modest, the amplitude of this deformation turns out to be comparable with the distortion of 
the elliptical surfaces of constant density at the CMB surface resulting from torsional oscillations 
in the core ( [Dumberry & Bloxham 20041 ). While this evidence alone does not imply a causal rela- 
tionship between the coseismic deformation field resulting from giant earthquakes and core flow 
perturbations, our results suggest that the CMB deformation of seismic origin has the potential 
to interfere with core dynamics. For instance, it has been recently a matter of debate whether a 
seismic perturbation of the CMB could trigger a flow instability, leading to a geomagnetic jerk 
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(Florind o et al. 2 005a; Dumberry 2005; Florind o~et al. 2005b~|) : in this case, a jerk should follow 
the seismic event, after a suitable time delay to allow for the signal to propagate through the 
weakly conducting mantle. 

Following the evidence that the main contribution to the deformation field comes from the 
lowest degrees coefficients of the spherical harmonic expansion, we have investigated the detailed 
time-dependent evolution of the perturbation to the elliptical part of the gravity field, J 2 . Our re- 
sults show that the Sumatra event produces a negative variation of J 2 , confirming the well-known 
tendency of earthquakes to reduce the Earth's oblateness (IChao & Gross 1987l|Alfonsi & Spada~9 98); 
this effect turns out to be further enhanced by the postseismic relaxation, according to our model. 

A sensitivity analysis, performed with varying asthenosphere viscosities, shows that, for values 
smaller than 10 18 Pa s, the postseismic effect on J 2 remains comparable with the main secular 
trend for several years after the event. In the next years the analysis of available J 2 data would 
allow us to put a lower bound to astheno spheric viscosity, which is still highly controversial issue 
(IPollitz et al. 1998UPiersanti 1999t|Marquart et al. 2005[). 



2 MODELING APPROACH 

To compute the coseismic and postseismic effects of the Sumatra earthquake on the core-mantle 
boundary, we used the semi-analytical model originally proposed by Piersanti et al. (119951) , which 
is a spherical, self-gravitating, incompressible model with Maxwell viscoelastic rheology. 

This model computes physical observables on the Earth's surface; however, its formulation 
straightforwardly allows to extract the deformation field at the CMB. In fact, the harmonic com- 
ponents of the physical quantities at the CMB (deformation, geopotential and stress tensor) are 
imposed as boundary conditions through the so-called "continuity matrix", so that the spheroidal 
and toroidal parts of the solution at the CMB (r = r c ) can be written as follows: 



y(r c ) 

Z ( r c) 



= I s (r c )c c 
= h(r c )c c 



4 V. Cannelli et al. 

where y is a 6-vector corresponding to the spheroidal part of the problem, z is a 2-vector 
corresponding to the toroidal part, I s and I t are the spheroidal and toroidal continuity matrices, 
whose expressions are given by Sabadini et al. (119821) and Piersanti et al. (1995), and c c and c c 
are, respectively, a vector and scalar constant to be determined by imposing traction-free boundary 
conditions at the Earth's surface. 

The perturbation to the gravity field elliptical term J 2 is related by definition to the I = 2, 
m = component of the geopotential as follows (Lambeck 1980): 

where Rr and M T are radius and mass of the Earth respectively, G is the gravitational constant 
and 2> o is the second-degree harmonic coefficient of the perturbation to the gravitational potential. 

While the computation of A J 2 involves the I = 2 harmonic term only, the evaluation of the 
deformation field requires the summation of hundreds of harmonic terms to gain a stable con- 
vergence. The CPU time needed to compute a single harmonic term increases strongly with the 
number of layers in the model (Bos chi et al. 20001) : in order to be able to employ a realistic, re- 
fined stratification and at the same time keep the computation time within reasonable limits, we 
adapted the analytical model formulation to the purely elastic and fluid cases by taking the limits 
t — > and t — > oo respectively, which in the Laplace domain correspond to s — > oo and s — > 0. 
In this way, we compute the full CMB deformation field in the elastic and fluid limits and give the 
transient postseismic evolution of A J 2 , which corresponds to Earth oblateness and is an indicator 
of the ellipticity evolution. 

The stratification model used in our computations is built by adopting the PREM (Preliminary 
Reference Earth Model) (Dziewonsk i & Anderson 1981b for mantle and crust and a uniform fluid 
core, with rigidity /i c = and density p c = 10.93 kg/m 3 , obtained by volume-averaging PREM 
core layers. The viscosity of the layers has been assigned by interpolating the viscosity model 
given by Mitrovica & Forte ( 2004 ). The resulting model has a total of 43 homogeneous layers and 
a uniform fluid core; its density, rigidity and viscosity profiles in mantle and crust are represented 
in figure [Q It is to note that, while the PREM model is compressible, the analytical formulation 



CMB and J2 perturbations from 2004 Sumatra earthquake 5 
of our model is based on an incompressible rheology; that is to say, we adopt a modified version 
of PREM with A — > 00. This approximation certainly affects our results, as discussed in detail 
by Nostra et al. (119991) , but presently it is an unavoidable choice if we want to take into account 
simultaneously viscoleasticity, self-gravitation and sphericity. 

The seismic source has been modeled using the five point sources obtained by Tsai et al. 
( 2005) . These seismic sources have been computed by fitting with the CMT method ( Dziewons ki et al. 1 981) 
the long-period seismograms from the IRIS Global Seismographic Network and account for a cu- 
mulative energy release corresponding to M w = 9.3. 

3 PERTURBATION OF THE CORE-MANTLE BOUNDARY 

In what follows we show and discuss the results we obtained in relation to CMB effects of the 
Sumatra earthquake. Figures [2] and [3] show, respectively, the scalar components of the dislocation 
vector u = (u r ,Ue, u<£) and the horizontal displacements at the CMB in an orthographic projec- 
tion, centered on the location of the composite CMT source obtained by Tsai et al. ( 2005 ), in the 
coseismic case. 

In particular, figure |2] shows the horizontal components along colatitude and longitude direc- 
tions (u0, u^), the radial (u r ) component and the absolute value of the displacement (|u|); the 
radial component shows a local CMB depression of about 4 mm by the event location. Starting at 
distances of 50° from the source location we observe a global vertical displacement of about 0.5 
mm. 

Figure [3] shows the horizontal displacements evaluated at CMB. The whole CMB surface is af- 
fected by appreciable displacements, with deformations still of the order of a fraction of millimeter 
even at extremely large epicentral distances. We also observe that the horizontal displacements are 
directed westward near the equator while near the poles the direction is opposite. The radial and 
horizontal displacement at the CMB computed by our numerical method are of the same order of 
magnitude. 

A spectral harmonic analysis has been performed in order to better understand the symmetric 
properties of the observed field deformations. Indeed, Dumberry & Bloxham (120041) pointed out 
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that only a CMB deformation satisfying axial and equatorial symmetry has the potential to interact 
with fluid core flows, possibly triggering a flow instability. 

Because of the location and north-south orientation of the fault plane, the CMB deformation 
field is approximately symmetric with respect to the equator plane (see figures |2] and |3]). On the 
other hand, the rupture geometry is not axisymmetric, so the deformation field has no intrinsic 
axial symmetry. Nevertheless, if we write the deformation field as a sum of spherical harmonic 
terms, we can extract the axial and equatorial symmetric terms and evaluate their amplitude. 

Let us now write the deformation field u(0, <p) as a sum of spherical harmonic functions: 

oo I 

u (^) = E E c i™yim{e,<i>) (i) 

1=0 m =-l 

with c\ m = (c[^, Cj^, c|^) being the vector whose elements are the harmonic coefficients of the 
expansion of the deformation along f, 9 and <p directions, respectively. The spherical harmonic 
functions Yj m are defined as: 



Yi m (e,d>) = J 2l + 1( f t m jv m ^ m (cosg) (2) 
47r (I + m)l 



with P im being the associated Legendre functions. The spherical harmonics satisfy axial and equa- 
torial symmetry only for even I and m = 0, so we can write the symmetric component of the 
deformation field as: 

u s (9) = cMO) (3) 

l even 

where we dropped the dependence on the m = spherical harmonic functions. The harmonic 
coefficients q , because of the orthonormality properties, can be immediately evaluated by: 

c, = J u(d,<p)Y l0 (6)dQ (4) 

with dQ = sm 9d9d(j) being the solid angle element. In figure H] we plot the harmonic coefficients 
QoforZ = 0,2, 20, computed by numerically integrating equation ©. The harmonic amplitudes 
show that a non-negligible amount of deformation associated with the lowest degrees satisfies the 
symmetry requirements. From figure |4] we see that the main contribution to symmetric term u s is 
given by the component of the CMB deformation. 

In figure [51 we show the symmetric term of the deformation u s , as in equation ©, and the 



CMB and J 2 perturbations from 2004 Sumatra earthquake 7 
associated residual = u — u s . Since the symmetric components of the deformation field 
are associated with lowest harmonic degrees, as shown in figure H they have been computed 
using only the spectral components / < 20. From figure [5] we see that the CMB deformation 
field exhibits spectral components that satisfy axial and equatorial symmetry accounting for a 
considerable part of the total deformation. For example, the component whose symmetric part has 
the largest amplitude is u^, which has a peak in the equatorial zone reaching 0.8 mm. Also the 
radial deformation, u r , shows a non-negligible symmetric component, with a considerable range 
of latitudes where the deformation amplitude exceeds 0.1 mm, while the component u e turns out 
to have the smallest symmetric term. 

We verified that the CMB deformation is strongly dependent from the dip angle of the seismic 
source: smaller dip angles result in less pronounced deformation effects on the CMB. This can be 
qualitatively explained by noting that in a source mechanism with a small dip angle there is little 
amount of slip in the radial direction, which is the component on which the CMB effects are most 
dependent. 

In figure [6] we plot the full CMB deformation field in the fluid limit (t — > 00). The horizontal 
components of the deformation field are greatly enhanced by the postseismic relaxation, while the 
radial component mean amplitude is comparable to the elastic case. As a result, the total deforma- 
tion vector u is about an order of magnitude greater than the elastic case, with peak values of a few 
centimeters. These results can be interpreted as a direct consequence of the viscoelastic relaxation 
of a Maxwell body, which in the fluid limit cannot sustain tangential stresses; therefore, the major 
effects are expected on the horizontal components of the deformation. Moreover, from a compar- 
ison of figures [2] and |6j we see that the deformation field in the fluid limit has a smoother spatial 
variation than the elastic limit, so we expect a further redistribution of the harmonic components 
towards lower wavelengths. 

4 COSEISMIC AND POSTSEISMIC EFFECTS ON J 2 

In this section we show the coseismic and postseismic effects of the Sumatra earthquake on the 
oblateness J 2 of the gravitational potential. This quantity is directly related to the Earth flattening 
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f = (a — c)/c, where c and a are the axial and equatorial radius respectively, so that positive 
variations of J 2 corresponds to an increase in the Earth's oblateness; for an homogeneous sphere, 
the simple relation J 2 = 2//5 holds. 

In figure[7]we show the long-term time dependence of J 2 resulting from the mass redistribution 
following the Sumatra earthquake. It was computed adopting the stratification model shown in fig- 
ure[TJ We see that, similarly to most thrust subduction earthquakes, the Sumatra event gives a neg- 



ative variation of J 2 , corresponding to a decrease of the Earth oblateness (Alfonsi & Spada 1998). 
This variation in the elastic limit is A J 2 = —0.30 x 10~ 10 and is further enhanced by the postseis- 
mic evolution, up to A J 2 = — 1.8 x 10~ 10 in the fluid limit. The transient evolution of A J 2 exhibits 
the largest variation for At ~ 10 2 4- 10 3 yr. This feature is to be ascribed to the detailed viscosity 
structure of the model, since it is dependent on the complex convolution of the model relaxation 
times. Incidentally, this kind of time dependence can be found in the postseismic modeling of other 
observables, such as the deformation field (INostro et al. 19991 IBoschi et al. 20 00). 

We note that our result for the elastic limit is in good agreement with that by Gross & Chao 
(|2006|) . which estimated AJ 2 = —0.24 x 10 -10 using the CMT source model obtained by Tsai 
et al. (2005) and a PREM elastic stratification. The coseismic variation of J 2 resulting from the 
Sumatra earthquake is therefore roughly equal to the variation occurring over a year due to the 
secular linear drift, which is j 2 ~ -0.28 x 10~ 10 yr" 1 ( [Cheng et al. 1997| ICox & Chao 20021) . 
and it is two orders of magnitude greater than the mean annual J 2 variation associated with global 
seismic activity (IChao & Gross 1987t|Alfonsi & Spada 1998) . 

This fact is particularly important since, as pointed out by Alfonsi & Spada (1998 ) the average 
effects of seismic activity and seismic tectonic movements tend to cancel each other both being of 
the order of 10~ 13 yr -1 and with opposite sign. Only with an exceptional event like the Sumatra 
earthquake, we have the chance to register its effects on J 2 . 

While the global coseismic deformation produces a jump in the J 2 evolution, the postseismic 
relaxation of the ductile astheno spheric layers is expected to give a continuous temporal variation 
of J 2 that will be superimposed to its secular drift. In what follows, we have computed the short- 
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timescale evolution of J 2 for various astheno sphere viscosities, to infer whether the viscoelastic 
relaxation may leave a detectable signature on the measured time-histories. 

For this purpose a simplified three-layer stratification model was employed, with an 80 km 
elastic lithosphere, a 200 km asthenosphere with variable viscosity and a uniform mantle with a 
constant viscosity of 10 21 Pa s. 

In figure [8] we show the time evolution of J 2 over a period of twenty years for asthenosphere 
viscosities 771 = 10 16 , i] 2 = 10 17 and i] 3 = 10 18 Pa s. We see that low asthenospheric viscosities 
yield very large variation rates in the first years after the event, as a result of the low associated 
Maxwell times = rjif/j,, with i = 1,2, 3. 

As such, these results can in principle be used to identify a lower limit for the asthenosphere 
viscosity on the basis of geodetic measurements of J 2 . In fact, if we assume a likely detectability 
threshold for deviations of J 2 from its secular drift and if no evidence of a such deviation is 
detected from available data, we can rule out the range of asthenosphere viscosities that produce 
perturbations on J 2 above that threshold. 

A reasonable value for the detectability threshold may be the associated formal error, which is 
about 10% of the measured value (Cheng et al. 1997 ), i.e. ~ 0.03 x lQ~ 10 yr~ l . In figure[8]the range 
of J 2 values below that threshold is represented by shaded area; as can be seen, for asthenospheric 
viscosities r] = 10 16 and 10 17 Pa s, the effect of the Sumatra event on J 2 would remain detectable 
for several years. 

In table \T\ the expected values of J 2 in 2005 and 2006 for asthenosphere viscosities ranging 
from 10 15 to 10 22 Pa s are given. As can be seen, from viscosities up to 10 17 Pa s we expect an 
evident signature in the data; viscosities greater than 10 19 Pa s should not produce a detectable 
signal, while the signal associated with a viscosity of the order of 10 18 Pa s lies marginal to the 
detectability threshold. We stress that this is just a general indication coming from a forward 
modeling; to apply this procedure to real data, a detailed sensitivity analysis will be needed. 
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5 CONCLUSIONS 

In this work we have shown that even the core-mantle boundary is affected by a significant amount 
of seismic deformation produced by the giant Sumatra earthquake, with coseismic radial displace- 
ments of the order of a fraction of millimeter over the whole CMB surface and horizontal displace- 
ments even larger. By analyzing the spectral components of the coseismic deformation field at the 
CMB surface, we found that most of this deformation is associated with low degree harmonics, 
and that the deformation field has considerable spectral components characterized by axial and 
equatorial symmetry. These symmetric components account for a radial deformation of the order 
of 0.1 mm and horizontal deformation with peak values slightly less than a millimeter. 

The CMB deformation field produced by the Sumatra earthquake turns out to be compara- 
ble to that resulting from core flow processes. In particular, we have verified that the CMB de- 
formation is of the same order of magnitude of that resulting from core torsional oscillations 
( [Dumberry & Bloxham 2004 ) and it is characterized by spectral components with similar symme- 



try. This suggests that the global deformation field from giant earthquakes has the potential to in- 
terfere with core processes. For instance, it has been recently suggested (Flori ndo et al. 2005a|) that 
a perturbation of the CMB topography of seismic origin could, at least in principle, trigger a core 
flow instability which can lead to a geomagnetic jerk. While this possibility is rather controversial 
( [Dumberry 2005 ; IFlorindo et al. 200551) . we have shown that the CMB deformation induced from 



a giant earthquake has a non-negligible amplitude and therefore should not be ruled out in playing 
a role in the triggering of core instabilities. In detail, Bloxham et al. (120021) argued that torsional os- 
cillations consistent with a geomagnetic jerk should have variations in amplitude of the order of 1 
km/yr. Our results show that the component of the CMB deformation axisymmetric and symmetric 
about the equator has an amplitude of the order of a fraction of millimeter. Although the amplitude 
of this deformation is modest, it is important to take into account that inside the core a torsional 
oscillations flow of 1 km/yr results from a distorsion of the elliptical surface of constant density 
of 0.2 mm at the Earth surface and of about 0.15 mm at the CMB ( [Dumberry & Bloxham 2004| ). 
This result opens the way to the sensibility that Sumatra event could really have triggered a jerk. 
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A definitive answer on this issue would require a comprehensive modeling of core-mantle interac- 
tion, that is beyond the scope of the present work. 

The postseismic evolution of the CMB deformation agrees with the known tendency of giant 
earthquakes to make the Earth rounder; that is, the net effect of global seismicity is a decrease 



of J 2 over time (Alfonsi & Spada 1998). In the fluid limit, only the horizontal components of the 
CMB deformation field are enhanced by about an order of magnitude, while the radial component 
remains of the same order of coseismic one; however, due to the extremely long time-scale of this 
processes, it is unlikely to expect any possible coupling with core flow geometry. 

An important feature of the long wavelength deformation field associated with this event is the 
sensitivity to astheno spheric viscosity. For asthenosphere viscosity values smaller than 10 18 Pa s, 
the postseismic J 2 variation lies above 10% of the mean secular trend for several years after the 
earthquake and therefore is expected to be detectable in the J 2 geodetic measurements. This would 
allow us, at least, to put a lower bound to mean astheno spheric viscosity that is still one of the 
outstanding issues in mantle rheology. In fact, when a detailed analysis of geodetic measurements 
of J 2 will become available, if a significant deviation of J 2 from its secular trend will be evident, it 
will give an indirect estimate of mean astheno spheric viscosity; on the other hand, if no deviation 
will be evident, it indicates a lower limit of rj > 10 18 Pa s. 
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Table 1. Expected J 2 annual variation rate during 2005 and 2006 as a function of asthenosphere viscosity. 
A simplified three layered stratification model was used. 



Asthenosphere viscosity 


j 2 2005 


j 2 2006 


(Pa s) 
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-0.003 


-0.003 


10 20 


-0.002 


-0.002 


10 21 


-0.002 


-0.002 


10 22 


-0.002 


-0.002 



Figure 1. Density, rigidity and viscosity profiles of the mantle and crust stratification model. The model 
assumes a uniform fluid core with density p c = 10.93 kg/m 3 , obtained by volume-averaging the corre- 
sponding PREM layers. 

P., Berger, J., Braitenberg, C, Van Camp, M., Lei, X., Sun, H., Xu, H., Rosat, S., 2005. Earth's Free 
Oscillations Excited by the 26 December 2004 Sumatra- Andaman Earthquake, Science, 308, 1139- 
1144. 

Piersanti, A., 1999. Postesismic deformation in Chile: Constrains on the asthenospheric viscosity, Geophys. 
Res. Lett, 26, 3157-3160. 

Piersanti, A., Spada, G., Sabadini, R., Bonafede, M., 1995. Global post-seismic deformation, Geophys. J. 
Int., 120, 544-566. 

Pollitz, EE, Burgmann, R., Romanowicz, B., 1998. Viscosity of oceanic asthenosphere inferred from remote 

triggering of eartquakes, Science, 280, 1245-1249. 
Sabadini, R., Yuen, D.A., Boschi, E., 1982. Polar wander and the forced responses of a rotating, multilay- 

ered, viscoelastic planet, Geophys. J. Int., 87, 2885-2903. 
Tsai, V.C., Nettles, M., Ekstrom, G., Dziewonski, A.M., 2005. Multiple CMT source analysis of the 2004 

Sumatra earthquake, Geophys. Res. Lett., 32, L17304. 



Figure 2. CMB coseismic displacement vector components along radial (u r ), colatitude (ug) and longitude 
(u^) directions and magnitude |u|. 



14 V. Cannelli et al. 



Figure 3. CMB coseismic horizontal displacements resulting from the Sumatra earthquake. The color scale 
represents the arrow length. 

Figure 4. Harmonic coefficients of the component of the CMB deformation field satisfying axial and equa- 
torial symmetry. 

Figure 5. Spectral component us of the CMB deformation field satisfying both axial and equatorial sym- 
metry and related residual ur = u — us- 

Figure 6. CMB displacement vector components in the fluid limit along radial (u r ), colatitude (ug) and 
longitude (u^) directions and magnitude |u|. 

Figure 7. Time-dependent postseismic evolution of the perturbation to the elliptical part of the gravitational 
potential, J 2- 

Figure 8. Time evolution of J2 variation rate for a three layered stratification model. Different lines rep- 
resent different values of asthenosphere viscosity. The shaded area corresponds to values of J2 below the 
detectability threshold (±0.03 x 10 _10 yr _1 ). 
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